Oscillation dynamics of multi-well condensates 
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We propose a new approach to the macroscopic dynamics of three- well Bose-Einstein condensates, 
giving particular emphasis to self-trapping and Josephson oscillations. Although these effects have 
been studied quite thoroughly in the mean-field approximation, a full quantum description is desir- 
able, since it avoids pathologies due to the nonlinear character of the mean-field equations. Using 
superpositions of quantum eigenstates, we construct various oscillation and trapping scenarios. 

PACS numbers: 03.75.Kk, 03.75.Lm, 03.65.Sq 



Interference phenomena belong to the most surprising 
features of quantum mechanics, which in many cases run 
against our classical everyday experience. Most spectac- 
ular is this nonclassical behavior, if it occurs in macro- 
scopic systems. Such cases can help to sharpen our eyes 
for quantum properties of matter. Therefore, it is always 
worth to point out situations, in which new examples of 
macroscopic nonclassical behavior might be accessible to 
laboratory experiments. The best candidates for macro- 
scopic quantum phenomena are supernuids and conden- 
sates. This is certainly one motivation for the enormous 
interest in condensates at the moment. In addition, con- 
densates serve as systems of analogy for many structures 
from very different subbranches of physics. This makes 
one hope that phenomena pointed out for condensates 
can find applications in various areas of physics. 

In this letter, we present a new idea for a quantum os- 
cillation of a condensate in three wells. It is a system, for 
which we have already presented a semiclassical analysis 
in [l[ . The complete classification of the eigenstates of the 
system in a dynamical representation makes it easy to se- 
lect appropriate eigenstates, which, when superimposed, 
result in a packet with the desired dynamical properties. 
We will show the following: By forming a packet from 
a few (in this letter we only use two) appropriate eigen- 
states, we produce oscillations between any pair of wells 
we want. The most counterintuitive case is an oscillation 
between the two outer wells, where the middle well is al- 
most empty and does not participate in the oscillations. 

Macroscopic effects like self-trapping and Josephson 
oscillations have been investigated for the two-well po- 
tential 0, 0, IH 0] and also for the three-well poten- 
tial [1, 13] • However, the behavior has been studied 
mostly in the mean-field approximation using the non- 
linear Gross-Pitaevskii equation or for specific quantum 
states. The subsequent investigation is done thoroughly 
in the many particle quantum system avoiding all the 
difficulties related to the nonlinearity of the mean-field 
equations. Since the described phenomena do not re- 
quire sophisticated preparation procedures, they should 
be directly accessible to experimental investigations like, 
e.g., in Refs. @,@]. 

We consider the situation, where a Bose-Einstein con- 



densate is confined in a linear configuration of three po- 
tential wells. In experiments, this can be achieved, e.g., 
by combining an optical dipole trap with a standing laser 
wave [8(. We restrict the analysis in this letter to three 
wells. A convenient set of basis functions for such sys- 
tems are the Wannier functions ip kt i(x) (loj . which are 
real-valued and exponentially localized in each well. The 
index k labels the potential sites, while I denotes the ex- 
citation in each well. Since we focus on condensed parti- 
cles, we can assume I = and skip the index in the fol- 
lowing. Thus, the three Wannier functions {tpi, tp 2 , ^3} 
span the one particle Hilbert space 7i\, For the sub- 
sequent discussion we only need the spacial localization 
properties of the basis, so that the analysis applies also 
to other physical situations like, e.g., atom chips 0. 

Expanding the field operator <3?(a;) in the one particle 
basis, &(x) — J2fk(x)ak, yields in the standard tight 
binding approximation the many particle Hamiltonian 
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The constants are given by 
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This is the celebrated Bose- Hubbard model [llj restricted 
to three sites. The coupling constant g describes the two- 
body interactions. In order to avoid degeneracies in e k , 
we assume a linear dependence on the additional param- 
eter A. In experimental realizations this can be achieved 
e.g. by a Stark field. The zero point of the energy can 
always be shifted in order to fulfill e = 0. We fix the 
values of these constants to A = 0.1, Ki 2 = K23 = 0.25 
and £ = 0.1 in order to describe the same situation as in 
[l[. These parameter values are experimentally accessi- 
ble with standard techniques and avoid degeneracies in 
the Hamiltonian ([I]) so that the present study shows the 
most general diversification. 
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We consider a system of N — 30 particles. The many 
particle Hilbert space has dimension L — (N + 1)(N + 
2)/2 and the eigenenergies E k and eigenstates k — 

1, L, can be calculated numerically exactly in the Fock 
basis |n) (h = (ni,n 2 ,n 3 )), 

I**) = E c ^> ( 5 ) 

with ^ = Cn^n2,n 3 £ R. In Q we introduced a 
dynamical (semiclassical) representation | Cp ) with 99 = 
(ipi, (f2, ^3) G T 3 in order to classify the eigenstates 
by the dynamical behavior of the corresponding classical 
system. We divided all eigenstates $k(0) into six cat- 
egories by assigning new geometrical quantum numbers 
to the states with respect to their localization properties 
on the toroidal configuration space T 3 . In the dynamical 
representation [12} , the eigenstates can be written as 

= <0l*k> = E (6) 

ni+ri2+n3— N 

Many eigenstates can be directly assigned to a category 
(A-El) with the idealized functional forms 

($\h,h)A = V 3 e tllVl e z{N - ll - l2)v2 e ll2tp3 , (7) 
{0% t) b = V 2 e^ 1 e^^+^s) Xr (^ 3 - H , (8) 
t) c = V 2 e i ^+ a ^ e a ^ X A<Pi - ¥> 2 ) , (9) 
t) d = V 2 e a ^ + ^' 2 e^ N - 1 ^ xA<Pi ~ <P»), (10) 
(<p\T d ,T a ) E1 = 2-V 2 V e iN ^+^y 3 

Xr d (<Pl + <P3 - 2^2) Xr a (<Pl - ^3) , (11) 

where rj = (2tt)~ 1 ^ 2 is a normalization constant, /, l\, I2, 
t, r a , Td € No are quantum numbers and Xn is a har- 
monic oscillator eigenfunction, which is localized in the 
configuration space [— it, 71"]. The sixth type (E2) con- 
sists of irregular superposition patterns and represents 
the chaotic regions of the classical phase space. We will 
not consider these states here, but study only the prop- 
erties of such eigenstates that are qualitatively similar to 
one of the functions of Eqs. ((7|)-(fTT|). In a good approxi- 
mation, the parameters on G No with ^ on = N — I only 
depend on the quantum number I. If an eigenstate sep- 
arates in the variable (fk, i.e. (fk enters only in the form 
e mfcVfc , then the condensate in well fc is decoupled from 
the rest and is an eigenstate of the number operator in 
semiclassical approximation (h — 1) 



I.e. it has a definite number of particles. The represen- 
tation (fl~2"j) of the number operators hk — a\a,k can be 
verified by comparing Eqs. (O and ([6]). So, hk acts as a 
momentum operator in the dynamical representation. 

In type (D), e.g., site 2 is decoupled from the rest and 
has the fixed number of particles N — I while sites 1 and 
3 are entangled with the total number of particles I. 



The overall picture is described as follows: There can- 
not be assigned global quantum numbers to all eigen- 
states due to the lack of symmetries in the system. How- 
ever, the whole Hilbert space can be divided into vari- 
ous sectors, where one of the sets of geometric quantum 
numbers of Eqs. ([?))- (fTTj) applies and characterizes the 
eigenstates completely. 
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FIG. 1: (Color online) Time evolution of the state |^} = 
(|0,0) B + |1,0) S )/V2. Shown are the particle numbers in 
well 1 (black, solid), in well 2 (blue, dashed), and in well 
3 (red, dash-dotted). The time is measured in multiples of 
T = 2tt/A. 

Now, we show how easy one can manipulate and under- 
stand the dynamical behavior of Bose-Einstein conden- 
sates using the dynamical representation (f7|)- (fTTj) . The 
easiest way to get a nontrivial dynamics is to superpose 
two eigenstates |$ a ) and of the system and study 
the time evolution: 

|*(t)> = c a e"^ 4 |$ ) + c b |$ 6 ) (13) 

with c a /b — |c a /b| e~ %la i b and a,b£ {1, ...,L}. The num- 
ber operators hk are fundamental observables for a Bose- 
Einstein condensate. Their time evolution is given by 
(n>k)t — (t)\ n k\^ (t)} , which can be calculated as 

(h k )t = \c a \ 2 (<S>a\hk\<S>a) + \c b \ 2 ($> b \h k \<S>b) 

+2|c Q | \cb\ |<$ |ft fc |$i>| cos(ftt + A 7 + £ fc ). (14) 

Here, = E a — E b , A 7 = ~f a — 7^ and we used the decom- 
position ($ |rifc|$&) = \(§a\hk\$b)\e l t k . The oscillation 
amplitude is proportional to the modulus of (& a \hk\$b), 
which we can calculate directly using Eqs. (|7])- (fTTj) . For 
the type (A) eigenstates, all three wells are decoupled, i.e. 
are eigenstates of all three number operators hk- There- 
fore, the non-diagonal matrix elements are zero and the 
behavior is time independent. Since the behavior of type 
(B) is equivalent to that of type (C), we discuss only 
types (C), (D) and (El). For convenience we choose 
\c a/b \ = 1/V2. 

Type (El) shows entanglement between all three wells. 
The matrix elements {r' a , Tl\hk\T a , 77,)^ can be calcu- 
lated analytically using Eqs. (Ilip and (fl"2"|) . First, the 
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integration region can be extended from the torus T 3 
to the the whole M 3 due to the localization of the os- 
cillator functions. Second, the matrix elements can 
be evaluated directly by transforming to new variables 
Xd = ipi - 2ip 2 + <p 3 , x a = cpi - tp 3 , x 3 = ip 2 , and rep- 
resenting the new differential operators —id/dx s by cre- 
ation and annihilation operators, 

- ^XrM = (b s - b\) X r s {x s ) (15) 

with s = d, a. The frequencies fld/a are given by the 
energy differences for the diagonal resp. antidiagonal de- 
gree of freedom of the local harmonic oscillators with en- 
ergies e^p, = Qd/a{] + 1/2)- The effective mass m c JJ a 
is the length scaling for the eigenfunctions (if\rd,T a ). 
The matrix elements of fik are different from zero only 




where the expectation value is calculated for the state 
\fy(t)). The dynamical representation of the eigenstates 
with the idealized wave functions (|T))- lfTTj) describes the 
dynamical behavior not only qualitatively but also quan- 
titatively. With m d s = 0.7 and £1^ = 1.8, the oscillation 
amplitude yields An^y 3 = 0.8, which coincides with the 
exact value obtained from Fig. [T] The different offsets 
of the oscillations in Fig. [1] are due to deviations of the 
phases of the exact eigenstates from the plane waves of 
the idealized wave functions (fTT|) . 

Beatings between eigenstates of type (C) are shown 
in Figs. [3 and [3] for states (|0,3) c + |0,4) c )/\/2 and 
(|24,4) c + |24,5) C )/V2 respectively. The matrix ele- 
ments are calculated along the same idea as in the pre- 
vious example: First, one transforms into new coordi- 
nates x\ — ip i — if2, X2 — </?3 — f2 and x 3 — if2, which 
maps the torus T 3 onto itself (shear). Then one extends 
the integration region in xi-direction to the whole K and 
represents —id/dxi as in Eq. (fT5|) . The matrix elements 
(I', r'\nk\l, t)c vanish for V ^ I and also for |r — t'\ > 1. 
So, the non-zero elements are (r = 1, 2, t' — t = 0, ±1) 

(l,r'\n r \l,T) c = a r ^ + (-1)^^-/)^/ — ^— ,(20) 
(l,T'\n 3 \l,T) c = l5 T ,, T . (21) 
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FIG. 2: (Color online) Time evolution of the state I*) = 
(|0,3)c + \0,4) C )/V2 and otherwise as in Fig.[T] 

if one quantum number coincides while the other differs 
by 0, ±1. We consider the case that r a — r' a . Then, 
the diagonal matrix elements with r' d = Td are given by 
(Td,T a \nk\Td,T a )Bi = N/3 for k = 1,2,3 and the non- 
diagonal elements yield for r' d = Td ± 1 

(T d ,T a \hi\Td,T a ) E i = isgn(r'-T) ^ dTd - , (16) 

(r d ,T a \n 2 \Td,T a ) E1 = -2isgn(r / -r) ]j d ^ , (17) 

i i \- \ \ r i \ I m< f f d mca 
(T d ,T a \n 3 \T d ,T a ) E1 = isgn(r -r) \J . (18) 

Here, f = max(r, r') and sgn(r) is the sign function with 
sgn(0) = 0. The particle numbers in wells 1 and 3 os- 
cillate in phase with the same amplitude, while ri2 oscil- 
lates with opposite phase. Fig. [1] shows an example for 
excitation numbers r d = 0, 1 and r a = 0. Shown is the 
numerically exact dynamics of the particle numbers nk 
defined by 

n k (t)= [ (&(x)$(x)) t dxn {a{d k ) t , (19) 

J well k 
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FIG. 3: (Color online) Time evolution of the state |^} = 
(|24, 4} c + 1 24, 5)o)/V2 and otherwise as in Fig. □ 

Fig. [2] shows particle oscillations between wells 1 and 2, 
while the number of particles in well 3 is constant with 
n 3 = 0.2. This is the effect of self-trapping directly shown 
in the quantum regime. The example shows the typical 
behavior of beatings in this class of eigenstates: A num- 
ber of n 3 = I particles is trapped in well 3, while the rest 
is oscillating coherently between the other wells with an 
amplitude proportional to the square root of the transver- 
sal quantum number. With m cff = 2.1 and Q = 1.8, the 
oscillation amplitude yields a value of An®y 2 = 2.7, which 
coincides with the exact value taken from Fig. O The 
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second example of this type in Fig. [3] shows the oppo- 
site regime where nearly all particles are localized in well 
3 and the rest oscillates with amplitude Arii/2 = 1-3. 
The matrix element yields for the amplitude the value 
AnfJ 2 = 1.2 (m oS = 1.8, ft = 0.3). 

In type (D), wells 1 and 3 are phase locked, although 
there does not exist a direct coupling term in the Hamil- 
tonian JT]). Both wells couple indirectly through the 
middle well, however also here we can completely de- 
couple the dynamics. The particle dynamics of the state 
(|4, 0)d + |4, 1) d)/V% is shown as an example in Fig. 2] 
In order to understand the behavior, we again have to 
evaluate the matrix elements (/', r'|nfe|/, r). The non- 
vanishing elements with V = I and |r' — r| < 1 can be 
written as (r = 1, 3, r' — r = 0, ±1) 

{I, T'\n r \l, r) D = t <5 r ,, T + z(2-r)sgn(r'-T)^^— ,(22) 

(l,T'\n 2 \l,T) D = (N -l)5 T ,, T . (23) 

Interesting are the two extremes, where all the popula- 
tion is in or outside of the middle well. In Fig.|4]nearly all 
particles (ri2 = 25.8) are concentrated in well 2 and the 
population of the others oscillates between them through 
the middle well. Calculation of the oscillation amplitudes 
by the matrix elements yields a value of An®? 3 = 0.46 
(m cS = 0.7, ft = 0.6) in comparison with the exact value 
of Ani/3 = 0.49. I.e., one particle is oscillating between 
the outer wells on average. The most impressive example 
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FIG. 5: (Color online) Time evolution of the state = 
(| 30, 2} D + |30,3)r,)/V2 and otherwise as in Fig. [TJ 



Concluding, we have shown various types of oscillation 
dynamics by superposing eigenstates. The classification 
of eigenstates in the dynamical representation provides 
a powerful tool in order to understand and manipulate 
the behavior of condensates in a few-well potential. For 
a larger number of particles N, there exist longer lad- 
ders of transversal excitations in each class leading to 
higher transversal quantum numbers. Accordingly, the 
amplitude of the particle oscillations can be increased by 
superimposing such higher excited states. There is also a 
generalization to longer chains of wells with oscillations 
between more distant wells, more on this in future work. 
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FIG. 4: (Color online) Time evolution of the state ]*} = 
(|4,0}d + |4, 1}d)/v / 2 and otherwise as in Fig.rjJ 

is shown in Fig. [5] The occupation of the middle well with 
ni = 0.2 is much smaller than the oscillation amplitude 
between the first and third well with An 1 / 3 — 1.3. Al- 
though the middle well is effectively not occupied, there 
is a particle oscillation of more than two particles through 
this well. Also here, the matrix element of the idealized 
wave functions yields a very accurate value for the oscil- 
lation amplitude of Anf /3 = 1.4 (m cff = 2.8, ft = 0.5). 
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